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ABSTRACT 

A  number  of  sequential  and  parallel  procedures  for  analyzing  the  response 
of  structural  complexes,  subjected  to  various  drives,  are  briefly  developed  and 
discussed.  The  impulse  response  vector  operator  is  defined  in  terms  of  impulse 
response  operators,  each  associated  with  a  unique  path  between  the  localized 
position  of  a  test  drive  and  a  localized  position  of  observation.  The  drive  is, 
correspondingly,  a  vector  and  the  response  is  the  scalar  product  of  the  impulse 
response  operator  and  the  drive  vectors.  A  sequential  procedure  of  subdividing  a 
structural  complex  into  a  number  of  coupled  dynamic  systems  is  stated  The 
formalism  is  then  stated  in  terms  of  matrices  and  vectors;  e.g.,  the  response  is  a 
vector,  each  element  represents  the  response  of  a  specific  dynamic  system,  the 
impulse  response  operator  is  a  matrix;  the  off-diagonal  elements  describe  the 
couplings  between  the  dynamic  systems,  etc.  If  the  dynamic  systems  are  chosen  so 
that  each,  in  isolation,  can  be  described  in  terms  of  an  eigen-impedance  operator, 
then,  in  addition,  a  modal  analysis  can  be  applied  to  the  multiple  dynamic  systems 
that  compose  the  model  of  the  structural  complex.  In  the  modal  analysis,  however, 
the  ranks  of  the  impulse  response  matrix,  the  response  vector,  and  the  drive  vector, 
are  swollen  by  the  modal  count,  usually  rendering  the  matrix  equation  for  the 
response  unwieldy.  The  modal  approach  may  be  substituted  by  a  wave  approach. 
In  this  parallel  approach,  the  propagations  in  the  dynamic  systems  are  described  by 
impulse  response  operators  that  are  commensurate  with  those  pertaining  to 
boundlessly  extrapolated  dynamic  systems.  The  finiteness  of  the  dynamic  systems 
are  accounted  for  by  junction  matrices;  a  junction  defines  the  boundaries  through 
which  dynamic  systems  interact  either  with  each  other  (transmissions)  or  with  self 
(reflections).  As  in  the  modal  approach,  in  this  wave  approach,  the  resulting 
formalism  is,  again,  rather  unwieldy.  It  is  shown  that  considerable  reductions  and 
simplifications  are  attained  if  the  complex  can  be  modeled  by  spatially  one¬ 
dimensional  dynamic  systems. 
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INTRODUCTION 


The  analysis  of  the  response  of  a  driven  complex  structure  is  relevant  not  only  to  dealing 
with  acoustical  problems,  but  also  to  dealing,  among  others,  with  optical,  biological,  and 
economical  problems.  This  generalization  and  encompassment  of  the  solutions  to  the  response  of 
structural  complexes  may  have  encouraged  wider  interest  and  contributed  to  the  variety  of 
approaches  to  the  analyses  of  these  generic  problems.  In  this  paper  the  acoustical  problems  are 
implicitly  epitomized. 

The  analysis  of  the  response  of  a  (structural)  complex  to  various  external  drives  is  often 
well  nigh  impossible.  One  must  then  resort  to  some  reduction  in  the  definitions  and  descriptions  of 
the  complex  and  the  external  drive  system  to  render  the  analysis  manageable.  A  reduction  is 
usually  achieved  by  modeling  the  complex  and  the  external  drive.  The  model  is  chosen  so  that  the 
analysis  can  be  performed.  Some  modeling  insists  on  details  so  that  one-to-one  correspondence 
between  the  response  of  the  model  and  the  actual  complex  is  substantially  preserved.  Measures 
and  devices  (active  and/or  passive)  that  are  instituted  to  control  the  response  (or  responses)  may  be 
compatibly  modeled.  The  elaborated  model  may  then  be  exercised  and  the  results  substantially 
related  to  the  controlled  complex.  In  such  applications  one  may  estimate  the  benefits  of  the 
instituted  controls  also  on  the  basis  of  one-to-one  correspondence.  On  the  other  hand,  some 
modeling  may  insist  merely  on  a  phenomenological  correspondence ;  one  may  be  satisfied  then 
merely  with  descriptions  of  the  response  that  manifest  a  phenomenon  (or  phenomena)  that  is 
common  to  the  simplified  model  and  the  actual  (structural)  complex.  Measures  and  devices 
(passive  and/or  active)  that  may  be  instituted  to  control  the  phenomenon  (or  phenomena)  may  then 
be  examined  on  the  elaborated  model  of  the  complex  and  phenomenologically  related,  by  analogy, 
to  the  controlled  complex.  In  such  applications  one  may  not,  however,  demand  precise  estimate  of 
the  benefits  accrued  from  the  institution  of  the  controls;  beneficial  trends  and  inclinations  is  all  that 
can  be  expected.  Moreover,  not  all  phenomenon  (or  phenomena)  can  be  accounted  for  in 
simplified  models  and,  therefore,  caution  should  be  exercised  in  any  phenomenological 
correspondence  procedures,  hi  this  paper  the  complex  is  variously  modeled  in  the  milieu  of 


phenomenological  correspondence.  Moreover,  it  is  not  modeled  with  the  intent  of  establishing  a 
particular  description,  but  rather  with  the  intent  of  exploring  analytical  techniques  and  procedures 
that  may  be  employed  in  deriving  general  phenomenological  descriptions.  Although  this  aim 
renders  die  paper  somewhat  abstract,  it  is  hoped  that,  nonetheless,  it  will  be  found  useful  by  some 
readers.  Finally,  the  relationship  is  established  between  modeling  that  utilizes  a  multiple  spatial 
dimensionality  of  the  actual  complexes  and  modeling  that  relies  solely  on  a  single  spatial 
dimensionality.  This  is  of  significance  to  the  modeling  used  in  recent  papers  [1-5]. 

The  basic  description  of  the  complex  is  cast  in  terms  of  an  impulse  response  function 
g(x  |  x',  1 1 1' )  so  that  the  external  drive  Pe  (x',  t'),  defined  at  the  position  { x',  t '},  will  generate  a 
response  p  (x,  t)  at  the  observation  position  {  x,  t);  i.e., 


I 


g(x|x',  t|t/)dx,dt,pe(x',t/)  =  p(x,  t)  ;  x=  £  sx,  ;  dx'^ndxj  , 

<-i  «-i  (la) 


where  s  is  a  unit  vector  lying  in  the  spatial  coordinate  x,  and  [s  n]  =  [n  s  ]  =  5^  -1  The  complex 
is  three-dimensional  when  s0  =  3,  two-dimensional  when  Sq  *  2,  and  one-dimensional  when 
Sq  =  1.  The  impulse  response  function  is  pure  if  it  is  a  functional  only  of  the  parameters  and 
quantities  that  describe  the  complex;  it  is  independent  of  the  drive  Pe  (x',  t')  and  the  response 
p  (x,  t ).  Situations  may  arise  in  which  it  may  be  convenient  and  instructive  to  express  equation 
(la)  in  spaces  other  than  {  x,  t};  e.g.,  in  the  {  x,  co  },{£,  t },  and  {£,  ©  }  -  spaces,  where  the 
wavevector  1c  and  the  frequency  co  are  the  Fourier  conjugates  of  the  spatial  vector  variable  x  and  the 


lThe  caret  is  used  herein  to  indicate  a  vectorial  representation.  A  caret  over  an  index  indicates  a 
" unit  vector";  e.g.,  £  in  equation  (la).  A  caret  over  a  variable,  e.g.,  x  in  equation  (la),  indicates 
that  it  may  be  multi-dimensional;  each  component  is  associated  with  a  given  normal  direction.  A 
caret  over  a  quantity,  e.g.,  g  in  equation  (8),  indicates  that  it  may  be  multi-path;  each  component  is 
associated  with  a  given  normal  path.  A  " unit  path  vector"  is  indicated  by  a  caret  over  the  index  of 
the  path;  e.g.  v  in  equation  (8). 


temporal  variable  t,  respectively.  The  spaces  of  convenience  may  be  chosen  on  die  basis  of  the 
externa]  drive,  the  response,  or  both.  For  example,  equation  (la)  may  be  expressed  in  the  forms 

J  g'(x|x',  t|o/)<^,d©/  Pc  (x'.cO  *  p  (x,  t)  ,  (lb) 

Jr(^i^®it/)^dt/pe(^,to-p(x,(B)  ,  (lc) 

JgCxIx'.mKo'ldx'dm'peCx'.wO'PC*.®)  .  ad) 


Jl/(^|k/,t|t/)dk/dt#fc(k/,t/)«p(tt)  , 


J  g'(k|fr,t|t')dfrdt'pB(fr,t/)-P(k.t)  , 


J  g(k|k',t|t')dk'dt'  pc(k',t')«p(£t)  ;  £-?sk,;  dk'-Ildk;  , 


jG(£|^m|©')d£dffl'Pe(lh<0«P(£«>)  . 


and  so  forth,  where  typically 


p(x,  a>)  -  ( 2ic)-l/2  J dt  p(£,  t)  exp  H<at)  , 


p(^,t)-(2Jcr(s°/2)Jdxp(x*0expa^x)  . 
|(x|$'»g>|oi>')-(2«)_1  JdtJ  dt'exp(-ia>t)g(x|£',  t|t')exp(i(DY)  . 


One  may  choose  equation  (lb)  if  the  external  drive  happens  to  be  highly  localized  in  the  frequency 
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domain  (steady  state  drive)  and  yet  the  response  is  desired  in  the  temporal  domain.  On  the  other 
hand,  if  the  external  drive  happens  to  be  highly  localized  in  the  temporal  domain  (transient  drive) 
and  one  desires  the  response  in  the  frequency  domain  one  may  choose  equation  (lc).  Further,  if 
the  external  drive  is  highly  localized  in  the  frequency  domain  and  the  response  is  desired  also  in 
this  domain,  equation  (Id)  may  be  chosen.  Other  combinations,  involving  also  the  spatial  and 
wavenumber  domains  can  be  similarly  classified.  However,  it  may  transpire  that  the  form  of  the 
impulse  response  function  may  dictate  in  part  which  of  the  spaces  are  to  be  conveniently  (and/or 
expeditiously,  within  a  Fourier  transformation,)  chosen.  For  example,  situations  may  arise  in 
which  the  impulse  response  function  may  be  stationary  in  one  and/or  the  other  dependent  variable; 
e.g.. 


g(x|x',  t|t')  =  (2jt)"1/2  g(x|x\  t-t')  , 


,-1/2  I 


g(iii',  tit') 

In  such  situations  one  finds  that 


g(x|x',  <o|G>')«g(x|x',  (0)  5 (©-co')  , 


I(£|£',t|0  =  I(£t|05(£-£'); 


So 


5(£-£')  =  n  8 &-£;)  , 


respectively,  where  typically 


g(x|x\  tt) 


(2rt)_1/2  Jdx  g  (x  |  x\  x)  exp (-icox) 


(3a) 

(3b) 


(3c) 


(3d) 


(3e) 


From  equations  (1)  and  (3)  it  is  apparent  that  stationarity  results  in  substantial  simplifications.  It 
reduces,  in  the  Fourier  transform  domain  that  corresponds  to  that  in  which  die  stationarity  is 
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present,  an  integral  operator  to  an  algebraic  factor.  The  availability  of  such  a  reduction  may  dictate 
the  choice  of  the  space  of  convenience  just  discussed.  Finally,  it  is  convenient,  for  die  most  part, 
to  deal  with  the  impulse  response  integral  operator  rather  than  with  the  impulse  response  function. 
The  impulse  response  operator  h(x  lx',  t  ( t')  is  related  to  the  impulse  response  function 
g(x|x',  tjt')  in  the  form 


e.g.. 


(4a) 


h(x|x',  t|t')pe(£',t')  a  J  g(x|£',  t|t')da'dt'pe($',t') 


(4b) 


where  pe  (x',t  0,  in  equation  (4b),  is  an  arbitrary,  but  well  behaved,  function  of  the  dependent 
vector  variable  (x' ,  t'}  Invariably,  except  for  a  few  simple  structural  complexes,  the  direct 
derivation  of  the  impulse  response  functions  or  operators  are  difficult  and  cumbersome.  One  must 
then  resort  to  modeling  procedures  that  will  give  the  derivation  a  chance.  In  this  paper  a  number  of 
analytical  approaches,  designed  to  derive  the  proper  impulse  response  functions  or  operators  of 
complex  structures,  are  considered. 

It  may  be  useful,  at  this  early  stage,  to  make  a  precis  of  die  story  of  die  paper,  hi  this  way 
some  idea  as  to  where  the  paper  goes  will  be  related  before  the  multitude  of  equations  and  the  maze 
of  notations  take  over.  The  response  of  a  (structural)  complex  that  is  excited  by  an  external  drive  is 
formulated  in  terms  of  an  impulse  response  function  and/or  operator.  The  operator  is  an  integral 
operator  that  spans  the  spatial  and  the  temporal  extents  of  the  complex  and  the  external  drive.  This 
formalism  is  depicted  in  equations  (1)  through  (4).  As  just  stated,  rarely  is  the  complex  and  the 
drive  sufficiently  simple  that  they  can  be  reasonably  represented  in  single  forms.  It  is,  therefore, 
proposed  to  model  the  complex  and  the  drive  so  that  the  description  can  be  decomposed  The 
decomposition  is  in  reference  to  unique  and  significant  paths  yielding  components  through  which 
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the  response  can  be  estimated  by  a  superposition.  The  hope  is  that  the  analysis  in  terms  of  these 
paths  may  lead  to  a  set  of  simpler  impulse  response  operators.  A  vectorial  arrangement  of  the  set, 
with  a  corresponding  set  of  drives,  may  render  the  formalism  simpler  in  certain  applications.  A 
more  drastic  step  in  the  analysis  is  then  taken.  The  complex  is  modeled  in  terms  of  elemental 
constituents  —  dynamic  systems,  each  of  which  is  amenable  to  a  reasonable  description.  The 
premise  is  that  the  impulse  response  operator  associated  with  each  dynamic  system  may  be 
reasonably  derived.  The  drive  is  similarly  modeled  in  terms  of  elemental  constituents,  each 
assigned  to  a  dynamic  system.  The  formalism,  as  expected,  is  cast  in  matrix  form:  The  complex  is 
described  in  terms  of  an  impulse  response  matrix  operator,  the  external  drive  in  terms  of  a 
compatible  vector,  and  the  response  is  then  a  vector,  each  element,  in  the  response  vector,  is  the 
response  of  a  specific  dynamic  system.  A  major  difficulty  in  this  analytical  modeling  is  the 
multiplicity  of  descriptions  required  and  the  requirement  for  the  description  of  the  interactions 
(couplings)  among  the  dynamic  systems.  Nonetheless,  this  modeling  technique  is  central  to  the 
paper.  This  modeling  technique  is  of  particular  appeal  in  those  situations  in  which  the  generic  and 
phenomenological  properties  of  the  response  of  a  collection  of  dynamic  systems  that  interact  may 
be  more  significant  than  the  precise  evaluation  of  that  response.  Nonetheless,  to  this  point  in  the 
development  of  the  analysis,  the  modeling  of  the  complex  in  terms  of  coupled  dynamic  systems  is 
formally  stated  and  the  resulting  equation  for  the  response  is  thus  of  limited  practical  utility.  The 
impotence  of  the  equation  is  evidenced,  for  example,  by  its  insensitivity  to  the  spatial 
dimensionality  of  the  dynamic  systems.  To  infuse  some  practicality  into  the  analysis  just 
developed,  two  major  approaches  are  taken:  1.  the  modal  approach,  2.  the  wave  approach. 

There  are  situations  in  which  the  modal  and  wave  approaches  can  be  readily  related.  The 
Sommerfield  Watson  transformation  constitute  such  a  relationship  [6,7].  However,  in  this  paper, 
the  relationship  is  not  a  paramount  topic.  Rather,  in  the  paper,  the  modal  and  the  wave  analytical 
procedures  are  considered  as  parallel,  but  distinct,  approaches  to  the  analysis  of  the  response  of 
models  that  cast  the  essential  characteristics  of  the  complex  in  terms  of  coupled  dynamic  systems, 
as  was  just  formally  recounted. 
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1.  In  the  modal  approach  the  couplings  are  expressed  in  terms  of  coupling  impedance 
operators  that  are  compatible  with  the  impedance  operators  of  the  individual  isolated  dynamic 
systems.  The  latter  impedance  operators  are  chosen  to  be  eigen-operators  so  that  each  defines  a  set 
of  eigenfunctions  and  a  corresponding  set  of  eigen-values.  The  sets  of  eigenfunctions  enable  one 
to  suppress  Ae  spatial  dependence  of  die  matrix  equation  for  the  response  vector,  thereby 
rendering  this  an  equation  of  temporal  dependence  only.  The  spatial  dependence  of  the  response 
can  be  recovered  if  the  modal  matrix  equation  for  the  temporal  portion  of  the  response  is  resolved 
and  the  eigenfunctions  are  explicitly  known.  It  emerges,  however,  that  the  modal  matrix  equation 
for  the  temporal  portion  of  the  response  is  of  higher  rank,  by  the  modal  count,  than  the  original 
equation.  It  appears,  therefore,  that  the  modal  approach  is  useful  when  the  modal  count  is  low  and 
the  eigenfunctions  are  Known.  Usually,  however,  the  increase  in  rank  and  the  general  difficulty  of 
specifying  the  eigenfunctions  for  most  dynamic  systems,  makes  the  modal  approach  impractical, 
unless  statistical  measures  are  brought  to  bear. 

2.  In  the  wave  approach  the  model  is  defined  in  terms  of  extrapolated  matrix  propagators, 
terminal  position  vectors,  and  junction  matrices.  A  propagator,  for  a  given  dynamic  system,  is 
defined  in  terms  of  propagation  in  which  back-scatterings,  from  regions  outside  the  insitu 
boundaries,  are  eliminated  by  an  appropriate  spatial  extrapolations  of  the  dynamic  system  The 
terminal  position  vectors  define  the  positions  of  the  boundaries  of  the  dynamic  systems.  The 
junction  matrices  define  the  boundary  conditions  and  the  couplings  among  the  dynamic  systems  at 
and  across  these  boundaries.  It  emerges,  from  the  derivation  of  the  matrix  equation  for  the 
response  in  this  format,  that  the  spatial  dimensionality  of  the  dynamic  systems  has  a  considerable 
influence  on  the  wave  formalism  The  interactions  via  the  various  paths  at  the  various  terminal 
positions  and  across  the  various  boundaries  lead  to  summations  (and  when  appropriate 
integrations)  galore;  so  much  so,  that  the  formalism  appears  too  hopelessly  interwoven  and 
entangled  to  be  processed  by  computers  of  reasonable  capacities,  today  a  useful  measure  for  the 
practir  lity  of  a  formalism  However,  dramatic  simplification  in  this  cumbersomeness  of  the 
formalism  results  when  the  dynamic  systems  are  one-dimensionaL  Further  meaningful 
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simplification  occurs  if  the  dynamic  systems  are  basic  —  the  propagation  in  the  dynamic  systems 
is  described  by  exronential-like  propagators.  The  response  of  those  complexes  that  may  be 
modeled  by  a  collection  of  coupled  basic  one-dimensicmal  dynamic  systems  is,  therefore, 
illustrated  and  is  now  readily  formulated  [1-5]. 

Finally,  it  is  proper  to  mention  the  comprehensive  review  by  Hodges  and  Woodhouse  on 
the  theories  of  noise  and  vibration  transmission  in  complex  structures  [8].  In  this  review  many 
aspects  of  the  subject  matter  covered  herein  are  discussed.  Indeed,  much  of  the  material  and 
references  presented  in  this  review  are  relevant  to  this  paper.  The  discussion  of  the  applicability  of 
the  formalism  of  the  acoustical  response  of  complex  structures  to  various  other  areas  of  physics,  is 
of  particular  interest  However,  the  intended  purposes  and  methodologies  of  the  review  do  not 
overlap  those  intended  and  developed  in  this  paper. 

MULTI-PATHS 

One  may  recognize  that  the  response  at  {x,  t>  may  be  contributed  to  by  a  number  of  possible  paths 
and  that  each  path  is  not  only  describable  by  an  independent  impulse  response  operator,  but  that 
only  a  portion  of  the  drive  may  participate  in  contributing  a  response  via  this  path.  The  complex  is 
then  modeled  by  a  number  of  paths  as  is  sketched  in  Figure  1.  With  the  aid  of  such  modeling  one 
may  state  that  for  each  path  the  equation  for  the  response  mimics  equation  (1)  so  that 

hv(x|x',  t|t')Pev(x',t')«Pv(x,  t)  ,  (5) 

where  v  designates  a  unique  path.  [cf.  equation  (4).]  It  is  helpful  and  instructive,  as  the  theoretical 
development  proceeds,  to  exemplify  and  contrast  the  derived  equations  with  those  pertaining  to  a 
simple  model.  Therefore,  the  derived  equations  will  be  followed  immediately  by  the 
corresponding  equations  that  are  appropriate  for  the  simple  model  of  the  complex.  The  simple 
model  may  also  be  employed  to  render  the  equations  more  explicit  Brief  references  to  this 
procedure  will  be  marie  at  the  appropriate  locations  in  the  text  Thus,  to  begin  this  procedure,  if  the 
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traversal  of  a  path  may  be  specified  merely  by  a  time  delay  tv,  the  impuse  response  operator 
hv  (x  |  x',  1 1 t')  in  equation  (5)  may  be  expressed  more  explicitly  in  the  form 


hv(x|x',  tlO  -  hv(x|x',  1 1  tO  A  Ctv)  , 


where  hv(x  |  x',  1 1 t')  plays  a  secondary  role  and  is  functionally  reasonable  and  available,  A  (t)  is 
a  temporal  operator  defined  by 

N  N 

A(t)f(t)«f(t-T)  ;  IlACXn)  -  A(SXn)  , 

n  a  (7a) 

the  path  v  is  assumed  unique,  the  subscript  and  superscript  (v)  identifies  the  path,  tv  designates  the 
elapsed  time  (retarded  time);  the  time  taken  by  a  response  induced  at  x'  to  transfer,  via  path  v,  and 
arrive  at  x,  and  finally,  f  is  an  arbitrary,  but  well  behaved,  function  of  time.  The  elapsed  time  t?  is 
a  function  of  x  and  x  and  is  independent  of  the  temporal  variables;  tv  «=  tv(x  |  x')  but  not 


tv  (x  |  x',  1 1 t')»  A  basic  example  of  the  impulse  response  operator  in  which  hv  (x  1  £',  1 1 t')  is 
more  explicitly  stated  is 


kv  r£  i  6'  ♦  i  ♦'> 


hv  (Xi x',  tit')  -  E*  (S | S',  t|t')  AW  ; 

F'  (x|x'f  t|t')  -  hv  (x|x')  Sfy-t')  ; 

hv(it|x/)-Jfv(x|x')d^  -  ;  5°*  -Jdft-t') dt'*< 

where  I*  (x  1  x')  is  independent  of  the  temporal  variables,  and 

At -f)  f(to-f(t) ;  n  -  AtN.i-to 


[cf.  equation  (4).]  To  cast  equation  (S)  in  the  format  of  equation  (1)  one  needs  to  construct  the 
impulse  response  operator  and  the  drive  in  vectorial  forms;  i.e., 

htfiS'.tio-r  h’'(4ii'.nt')v  , 
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$,(«',  t')  =  2 t')  ;  [0S]  =  [&«]  =  5,„  , 

O 


(9) 


where  again  the  caret  atop  a  quantity  signifies  its  vectorial  nature.1>2  [cf.  equations  (la)  and  (lg).) 
For  die  simple  model,  equation  (8a)  assumes  the  form 


h(x|x',  tlO-Ih^xlx',  t|t')  A(tv)v  . 


(8b) 


The  grand  equation  for  the  response  is  then  derived  from  equations  (1).  (4),  (8),  and  (9)  in  the 
form 


[h($|^t|Ofc(£',t')]=p(a,t)  ;  p&O-Spv&t)  . 


(10) 


It  is  noted  that  a  path  is  defined  either  in  terms  of  components  in  the  impulse  response  operator,  in 
the  drive,  or  in  both.  Although  not  exactly  in  this  format,  equations  (5)  and  (10)  have  proved 
useful  in  handling  problems  relating  to  structural  response.  Problems  relating  to  propagations  at 
sea  have  made  greater  use  of  such  decompositions  than  structural  vibrations  have.  The  reason  for 
the  limited  deployment  of  this  formalism  in  solving  problems  related  to  structural  vibrations,  stems 
from  the  difficulty  and  cumbersomeness  that  often  still  remains  in  the  description  of  the  partial 
impulse  response  operator  hv(x  |  x',  1 1 t')  — the  impulse  response  operator  describing  the  transfer 
via  a  single  unique  path;  path  v.  This  statement  holds  true  even  for  die  partial  impulse  response 
operator  of  models  that  are  simple.  [The  simplicity  of  these  models  may  arise  because  they 
represent  simple  practical  structural  complexes  and  not  necessarily  because  they  represent 
excessively  simplified  complexes.]  Much  of  the  progress  made  has  been  linked  with  statistical 
measures  designed  to  perform  the  summation  over  die  paths;  e.g.,  tty  ensemble  averaging.  In 


2Paths  (and,  when  approapriate,  positions  at  junctions)  are  cast  in  terms  of  sums  over  discrete 
values.  However,  in  principle  these  sums  may  be  continuous  so  that  the  summations  may  be 
appropriately  replaced,  in  part  or  in  whole,  by  integrations.  This  transition  from  the  discrete  to  the 
continuous  may  be  interpreted  from  ray,  via  beam,  to  wave  descriptions.  In  this  sense,  in  the  text, 
the  reference  to  wave  propagation  is  more  akin  with  ray  propagation. 
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these  cases  die  partial  impulse  response  operators  need  not  be  precisely  defined  and  described. 
Averaging  of  this  kind  has  been  applied  with  some  success  to  limited  and  simplified  situations  in 
structural  vibrations;  e.g.  References  9  and  10.  On  the  other  hand,  averaging  of  this  kind  has  been 
traditionally  quite  acceptable  in  solving  problems  relating  to  propagations  at  sea;  e.g.  References  11 
and  12.  On  a  typical  basis,  statistical  measures  may  yield  from  equation  (10)  a  simplified 
equation  (1).  However,  in  this  paper,  rather  than  investigate  die  restrictions  that  are  imposed  on, 
and  explore  the  advantages  that  are  accrued  from  the  resulting  equation,  one  may  inquire  whether 
alternate  and/or  further  analytical  decompositions  exist  from  which  more  profound  benefits  may  be 
attained. 


MULTI-DYNAMIC  SYSTEMS 

One  may  attempt  to  subdivide  the  complex  into  a  number  of  coupled  dynamic  systems  as  is 
sketched  in  Figure  2.  [cf.  Figure  1.]  A  contrived  model  may  then  be  constructed  in  which  the 
dynamic  systems  are  rendered  uncoupled.  The  uncoupling  of  the  dynamic  systems  may  be 
instituted,  for  example,  by  artificially  isolating  or  extrapdating  the  dynamic  systems.  In  this 
contrived  model  each  dynamic  system  admits  to  an  impulse  response  function  that  is  hopefully 
reasonable  to  describe.  The  equation  for  the  response  of  the  contrived  model  of  the  complex,  in 
which  the  dynamic  systems  are  rendered,  artificially  if  necessary,  uncoupled,  may  be  cast  in  the 
form 


h®  (£l*'> 1 I1')  P«  (*'»  O  -  Pu  (£»  *) 


(11a) 


where 


so 

I 

i 


s  xsj 


dx; 


S0  A 

n<^ 


(12) 
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h„($|$',t|t')  -  (hujC^jl^.tlOSji) 


(13a) 


Pe(x'>  t')  «  {pci(x'i ,  t')}  ;  pu(x  ,  t)  -  {p^j(Xj ,  t)}  . 

-  ~  ~  (14a) 


P«  ( xj »  O  is  the  drive  applied  at  positions  {xj,  t' }  to  the  (i)th  dynamic  system, 
hnj  (xj  |  Xj ,  1 1 1')  is  the  impulse  response  operator  of  the  uncoupled  (j)th  dynamic  system,  and 
Puj  (xj ,  t)  is  the  response  at  position  {xj ,  t }  in  the  (j)th  dynamic  system.  In  component  form 
equation  (11a)  is 


hnj(Xj  |  Xj,  t  1 1  )  Pej(  Xj  ,  t  )  *  Puj(Xj  ,  t  ) 


(15a) 


It  is  noted  that  die  dynamic  systems  need  not  span  the  same  spatial  dimensionality.  Thus,  Xj  may 
be  of  different  spatial  dimensionality  than  that  of  jq  with  i*j. 


The  insitu  model  of  the  complex  must  account  for  the  various  couplings  that  may  exist 
between  the  multiple  dynamic  systems.  Indeed,  the  specification  of  the  couplings  is  crucial  in  this 
modeling.  It  is  proposed  that,  in  turn,  the  modeling  may  reveal  the  manner  in  which  the  couplings 
need  be  specified.  This  interaction  between  the  modeling  and  die  specification  of  the  couplings  is 
essential  to  the  success  of  this  formalism.  [It  is  to  be  understood  that  in  the  specification  of  die 
couplings  are  included  deviations  in  the  boundary  conditions  under  which  equation  (1  la)  is  stated. 
Therefore,  the  specification  of  the  couplings  may  cause  changes  in  the  definitions  of  some  dynamic 
systems  without,  necessarily,  causing  changes  in  the  fact  that  they  are  uncoupled;  Le.,  some 
dynamic  systems  may  be  uncoupled  in  the  insitu  model]  Taking  account  of  the  couplings,  one 
may  formally  state  die  equation  for  the  response  in  die  matrix  form 

h(i|j*',tr)p.(i'.t')-p(i,t)  ,  (lib) 


or  equivalently  in  the  component  form 
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?  hji ( xj  |  xj,  t  It') Pci  (  xj ,  t')»  P|j  (xj ,  t )  , 


(15b) 


where 

h(x|x\t|t')  «  (hji^jU'i.tlO)  , 

Pe(x t')  »  {p*  (xj  ,  t')}  ;  P  (£  >  t)  «  {«  ,  0}  . 

mt  0t*  ** 


(13b) 

(14b) 


[cf.  equations  (11a),  (15a),  (13a),  and  (14a),  respectively.]  In  these  equations,  the  impulse 
response  operator  hji  (Xj  |  xj ,  1 1 1')  is  the  "transfer  operator"  from  the  drive  position  at  {xj ,  t' } 
in  the  (i)th  dynamic  system  to  the  observation  position  at  {  Xj ,  t }  in  the  (j)th  dynamic  system, 
Pei  ( xj ,  t')  is  the  external  drive  at  {xj ,  t' }  on  the  (i)th  dynamic  system,  and  Pj  (Xj ,  t)  is  the 
response  at  {xj ,  t }  in  the  (j)th  dynamic  system  due  to  all  the  external  drives  in  the  various 
dynamic  systems.  One  may  also  cast  equation  (1  lb)  in  the  more  decomposed  form 


(11c) 


where  the  (external)  drive  matrix  and  the  response  matrix  are  defined 

Pe(x',t')»(pei(x[,  t')8ik)  ;  P(x  ,t)  -  (pji(Xj.t)) 

m  1  m 


(14c) 


respectively.  In  component  form  equation  (1  lc)  reads 

hji($j  I  xj,  1 1 0  pei  (  xj ,  t')  -  pu  (xj ,  t )  .  (15c) 


[cf.  equations  (15a)  and  (15b).]  It  is  obvious  that  when  the  dynamic  systems  are  artificially 
rendered  more  and  more  uncoupled,  in  the  prescribed  manner,  the  result  is 


m  ~  ~  m  ~ 


(16) 
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The  inverse  of  equation  (16)  is  mote  difficult  and  cumbersome  to  establish.  Not  only  off  diagonal 
elements  appear,  rather  than  disappear,  but  the  diagonal  elements  undergo  changes.  Therefore,  the 
establishment  of  the  uncoupled  impulse  response  matrix  operator  htt  ( x  1  x',  1 1 1' )  may  not  imply  a 
direct  and  an  immediate  step  in  the  process  of  deriving  the  impulse  response  matrix  operator 
h  (x  |  x',  1 1 1' ).  Nonetheless,  situations  may  arise  in  which  such  an  establishment  may  constitute  a 
t  first  step  in  the  process  of  deriving  h;  e.g.,  hu  may  be  considered  the  zeroth  order  approximation  to 

h.  In  this  context  modeling  the  complex  in  terms  of  multiple  dynamic  systems  necessitates 
specifying  the  couplings  among  the  dynamic  systems  in  a  manner  that  is  compatible  with  the 
definition  of  the  contrived  model  of  uncoupled  dynamic  systems. 

The  derivation  of  the  impulse  response  matrix  operator  h  (x  |  x 1 1 1' )  for  the  subdivided 
complex  can  be  achieved  by  a  number  of  different  approaches.  Invariably  in  these  approaches 
further  analytical  decompositions  are  instituted.  Are  these  further  analytical  decompositions  of  help 
,  in  reducing  the  difficulty  and  cumbersomeness  of  ascertaining  the  response  of  structural  complexes 

excited  by  complex  external  driving  systems?  Two  distinct  approaches  are  subsequently 
discussed,  and  comments  are  made  with  respect  to  the  answer  to  this  question. 


MODAL  APPROACH 

An  analytical  decomposition  that  is  widely  familiar  is  the  modal  analysis.  A  model  in 
which  the  multiple  dynamic  systems  are  rendered  uncoupled  by  isolation  is  amenable,  by  its  very 
structure,  to  modal  analysis.  Indeed,  the  isolated  dynamic  systems  in  the  complex  are  conveniently 
chosen  in  a  manner  that  make  them,  individually  and  collectively,  suitable  to  a  modal  analysis. 

The  model  is  sketched  in  Figure  3.  [cf.  Figure  1  and  2.]  The  equation  for  the  response  is  then 
typified  by  equation  (1  la).  It  is  assumed  that  each  appropriately  isolated  dynamic  system  is 
characterized  by  an  impedance  operator.  The  impedance  operator  Zj  (xj,  t)  is  the  inverse  of  the 
impulse  response  operator  huj($j  (  $j,  t  |t')of  the  (j)th  dynamic  system,  [cf.  equation  (13a).]  It 
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is  further  assumed  that  the  chosen  impedance  operators  are  eigen-impedance  operators.  It  is  noted 

that  the  definition  of  the  eigen-impedance  operators  is  not  unique;  there  are  usually  a  number  of 

ways  in  which  the  dynamic  systems  in  die  complex  may  be  isolated.  A  convenient  set  of  eigen- 

impendance  operators  is  then  selected  to  define  the  model  of  the  complex.  Each  eigen-impedance 

operator  defines  a  closed  (complete)  set  of  orthogonal  (normal)  eigenfunctions  and  a  corresponding 

set  of  eigen-values.  The  modal  (eigen-)  equations  for  die  response,  the  closure,  and  orthogorality  | 

are  expressed  by 


Zj(Xj,  t)  <t>j„(Xj)  -  Zja(t)  4>jn(Xj)  , 

(17a) 

?  $jn(^j )  <M*j  )  -  8(Xj  -  xj)  ; 

(17b) 

1  $jn(Xj)  dXj  <(>jm(Xj)  *=  Sum  » 

(17c) 

respectively,  where  Zj  (xj,  t)  designates  the  eigen-impedance  operator,  <t>jn(£j )  designates  an 
eigenfunction,  and  z^ft)  designates  an  eigen-value  for  the  (j)th  dynamic  system,  n  (and  m)  is  die 
modal  (eigen-)  index,  and 

5<xj - i; )  -  ns<x^- .  (18) 

Strictly  the  modal  index  need  be  cast  in  terms  of  a  vector.  That  is,  die  index  n  (and  m)  in  die 

a  s°  A 

equations  in  this  section  need  be  replaced  by  die  vector  index  n  (and  m),  where  n  •  5^  sn,.  Since 

the  modal  approach  is  widely  familiar,  such,  and  other,  details  in  the  modal  approach  are  not 
dwelled  upon  in  this  paper.  From  equations  (11),  (14),  and  (17)  one  may  state  and  readily  obtain 

pri(x'i,  t')  *  5  Pwm(t')  WxJ)  ; 

Pj(xj,t)  *  $  Pjn(t)  <M$j)  ’  (19a) 

where 
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Peim(t')  *  J  Pei(xj,  t')  dxj  ; 

Pjn(t)  -  J  Pjn(Xj.t)  dXj  <|>jn(Xj )  , 


(19b) 


The  insitu  model  is  now  constructed  by  introducing  the  couplings  among  die  dynamic  systems. 
The  couplings  are  compatibly  defined  in  terms  of  coupling  impedances.  An  impedance  matrix 
equation  for  die  response  may  then  be  defined  in  die  form 

zp(x,  t)«  pe(x,t)  ;  p(x,  t)  *  {pj(xj,  t)}  ;  pe($.  0  *  Pejdtj.  t)  ; 

i»(jj(Sj, t)Sji-| 2ji(ijix;,t)  d£;  J SfXj-xi)  d4i—)  , 

(20) 

where  Zj  (xj,  t )  is  as  defined  in  equation  (17)  and  Zji  (*j  I  xj ,  t )  is  the  coupling  impedance,  this 
impedance  describes  the  coupling  between  the  (i)th  and  the  0th  dynamic  systems.  [It  is  noted  that 
to  render  the  dynamic  systems  isolated;  equation  (16),  one  merely  removes  the  coupling 
impedances;  or  more  generally,  one  renders  z  diagonal.]  Transforming  equation  (20)  into  the 
modal  domain  and  making  use  of  die  orthogonality  of  the  eigenfunctions  yields 


z(t)  £(t)  «  pc(t)  , 


(21a) 


Pe(0  *  {•••Pei(t,)”*3’  J  PeiO*)  *  {pein^O}  i 
p(t)  »  {•••£j(t)---}  ;  £j(t)  «  {Pjn(t)}  . 

Z  *  (zjn(t)  5ji  5nm  —  ZjnimO))  * 

where  z^,(t)  is  defined  in  equation  (17a)  and 

Zjnia(t)  «  J  $jn(*j)  d^j  *ji(*jl*i.O  dx|  W*!)  • 


(22a) 

(22b) 


(23a) 


or  equivalently 


Zj| (Xj ) Xj,  t)  *  jf-  $jn(Xj)  zjnim(t)  $im(Xj ) 


(23b) 


The  quantity  2^(0  is  the  coupling  eigen-value  describing  the  modal  coupling  between  the  (m)th 
mode  in  the  (i)th  dynamic  system  and  the  (n)th  mode  in  the  (j)th  dynamic  system.  In  this 
description,  is  the  self-coupling  eigen-value;  it  is  these  eigen-values  that  modify  the  diagonal 

elements  when  couplings  are  introduced.  Equation  (21)  may  be  inverted  so  that  the  equation  for 
the  response  is  expressed  in  terms  of  the  impulse  response  matrix  operator  h  that  is  associated  with 

z.  The  inverted  equation,  in  matrix  form,  is 

P(t)«h(t|t')  pe(t')  ;  h(t|0  =  J  g(t|t')dt'  ••••  ; 

h(t|t  )«  (hjoiie ( 1 1 1  ))  ;  g ( 1 1 1  )*  (gjnim(t|t  ))  ,  (21b) 


or  equivalently,  in  component  form,  is 


5  2  hjnini(t  1 1 /)  pens(t  )  m  pja(0  » 
hjnijn(t  |t')  *i  gjnimft  !t/)dt/"»*  , 


(21c) 


where 


gjnim(t  1 t')  «  J  <|>jn(Xj)  dxj  gji(Xj  |  x[,  1 1  t')  dxj  tem(xj) 


(24a) 


or  equivalently 


gji (Xj  |  Xj,  1 1  t')  *  ^  J  ^jnCXj )  gjnim(t  1 1 ')  ^hmfXj ) 


(24b) 


(cf.  equations  (15b)  and  (23).]  The  inversion  of  equation  (21a)  into  equation  (21b),  and  to  that 
matter,  vice  versa,  is  greatly  facilitated  if  the  dynamic  systems  and  die  couplings  among  them 
remain  unchanged  with  time  —  the  dynamic  systems  and  the  couplings  among  them  are 
temporarily  stationary.  That  is,  z(t)  is  pure  differential  operator  in  the  temporal  domain  and 


g  (t  1 1')  »  (2it)*1/2  g  (t  - 1").  [cf.  equation  (3a).]  Assuming  this  temporal  stationarity  and 
performing  Fourier  transformations  with  respect  to  die  temporal  variable  on  equations  (21a), 
(21b),  (23),  and  (24)  yield 


|«D)  g«D)«ge(CD)  ;  I(a»«(2jttim(G»)  , 


(25a) 


ff(o»-|(u>)ge(fi»  ;  |(a>) *  [|(©)]'1 


(25b) 


zji(xj  I  f  §  9jn(*j)  <i>im(*i) 


(23c) 


8ji(^j  I  ®)  *  £  §  4>jn(*j)  IjBim(0>)  fom(*i) 


(24c) 


respectively,  where  ff(  a>),  §e(  G>),  and  f  ( co)  are  the  Fourier  transforms  of  p  ( 0,  pe(  t),  and 
(2jc ym  g(t-t') ,  respectively,  and  00)  is  the  Fourier  eigen-value  matrix  of  die  temporal  matrix 
operator  z(t); namely, 

z  ( t)  (ate)*112  exp  (icot)  5ji  5nm)  -  2(»)  ((2* )- 1/2  exp  (icot)  %  Snm)  .  „ 


It  is  noted  that  die  inversion  from  equation  (25a)  into  equation  (25b),  and  vice  versa,  is  algebraic. 
The  simplicity  of  the  inversion,  in  the  frequency  domain,  when  a  temporal  stationarity  prevails,  is 
thus  clearly  demonstrated.  However,  as  equation  (25a)  indicates,  die  formalism  just  discussed  is 
predicated  not  only  on  die  intimate  knowledge  of  the  impedance  operators  that  describe  the  isolated 
dynamic  systems,  but  also  on  knowledge  of  the  coupling  impedance  operators.  Such  knowledge 
is  rarefy  available.  In  this  connection  it  may  be  noted  that  the  eigenfunctions  need  not  be  directly 
defined  by  an  actual  impedance  operator.  For  example,  eigenfunctions  that  are  geometrically 
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defined  are  quite  acceptable  in  deriving  the  equations  in  this  section.  However,  the  geometry,  like 
the  spatial  dependence  of  the  eigen-impedance  operators,  needs  to  remain  stationary  with  respect  to 
the  temporal  variable.  Using  artificial  eigenfunctions  is  not  without  penalty,  it  renders  the  modal 
quantities  and  parameters  also  artificial  and,  therefore,  not  physically  interpretable. 

Comparing  equation  (15)  with  equation  (21)  it  is  observed  that  the  spatial  dependence  is 
stripped  off  by  the  modal  procedure;  a  simplification  indeed,  however,  at  some  expense.  The  rank 
of  the  matrix  h  is  swelled  from  that  of  h ,  defined  in  equation  (13),  by  the  number  of  modes  within 

— *  at 

the  complex.  Even  in  simple  complexes,  such  swelling  may  be  quite  substantial  and  equation  (21) 
and  (25)  may,  thus,  become  unwieldy.  In  addition,  seldom  is  knowledge  of  the  natural 
eigenfunctions  readily  available,  even  in  these  simple  complexes.  Means  to  reduce  the  rank  of 
equation  (21)  and  (25)  may  prove  essential  and  the  requirement  for  knowledge  of  the 
eigenfunctions  (and  the  eigen-values  of  the  coupling  impedance  operators)  may  need  relaxation. 
Again,  statistical  means  may  have  to  be  devised  to  induce  such  a  reduction  and  to  suppress  the 
need  to  know  detail  forms  of  the  eigenfunctions  of  the  isolated  dynamic  systems  and  the  eigen¬ 
values  that  specify  the  couplings  among  the  dynamic  systems  of  the  complex.  The  statistical 
energy  analysis  (SEA)  is  one  such  a  device  [13, 14].  The  rendering  of  SEA  from  equation  (21) 
and/or  (25)  lies,  however,  beyond  the  scope  of  this  paper. 


WAVE  APPROACH 

Another  approach  to  an  analytical  decomposition  is  a  wave  propagation  approach.  This 
approach  is  net  new;  it  has  been  employed  in  die  vibrational  analysis  of  complex  structures  [8]. 
Indeed,  in  part,  this  approach  is  employed  in  Section  H  In  this  paper,  however,  an  attempt  is 
made  to  systematize  this  approach  and  finally  recover  and  extend  some  recent  work  by  the  authors 
[4,5].  Equations  (5)  and  (1 1)  are  considered  basic.  In  this  approach  die  dynamic  systems  in  the 
contrived  model  are  rendered  uncoupled  by  extrapolation.  An  uncoupled  dynamic  system  in  the 
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complex  is  chosen  so  that  a  propagation  mode  can  be  assigned  to  it;  namely,  were  the  (j)th  dynamic 
system  extrapolated  boundlessly,  the  propagating  path  v  in  it  can  be  defined  in  terms  of  the 
propagator  function  h^j  ( xj  |  xj ,  t  1 1')  so  that  a  corresponding  external  drive  Pevj  (xj,  t')  would 
induce  a  response  P«j  ( Xj ,  t )  at  the  position  {  xj ,  t  }  on  the  extrapolated  dynamic  system .  The 
analytical  expression  for  this  process  is 

h^(xj lx/,  t| t')  pevj(^; ,  to  =  p^(xj ,  t)  . 


[cf.  equations  (5)  and  (15a).]  Again,  were  the  path  of  propagation  a  simple  path,  in  equation  (27) 
hXj(xj  |  Xj  ,  1 1  tO  assumes  the  more  explicit  form 


hXj(^j  |  x/  ,  tltO-h^tfjlSj',  t|t')Aj  , 


rv  /A  .  a  * 


Aj  s  A(x£) 


(28) 

(29) 


where  the  quantity  x-  designates  the  time  taken  by  a  response  induced  at  xj  to  transfer,  via  path  v 
on  the  extrapolated  dynamic  system,  and  arrive  at  Xj ;  (x3 1  Xj') .  [cf.  equation  (5) 

through  (7).]  The  extrapolation  is  considered  proper  provided  it  is  instituted  in  a  manner  that 


would  avoid  back-scattering,  from  regions  outside,  into  regions  occupied  by  the  actual  dynamic 
system.  It  is  noted  that  the  extrapolation  is  not  unique;  there  are  several  ways  in  which  it  can  be 
instituted.  However,  an  extrapolation  of  this  kind  is  essential  and  central  to  the  wave  approach 
here  proposed.  Equation  (27)  for  the  multitude  of  dynamic  systems  may  be  stated  in  the  matrix 
form 


where 


i'.tIO  -  pl(x,t) 


(30a) 
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.<£l£'.*lO  «  (h^j  (Xj  I  xj,  t  I  t')  Sji) 


(31a) 


Pcv(x',0-{pevj(x;,t')}  ;  pV(£,t)  «  {p^tfj.t)}  . 

-  (32) 

In  the  wave  approach,  equation  (30)  is  equivalent  to  that  of  the  generic  expression  specified  in 
equation  (1  la).  [The  prescribed  "extrapolation"  in  the  wave  approach  plays  the  role  of  die 
"isolation"  in  the  modal  approach  discussed  in  the  preceding  section.  These  two  prescriptions  are 
designed  to  derive  equation  (1  la)  in  forms  that  are  fundamental  to  each  of  the  two  approaches.] 
From  equations  (27),  (28),  and  (32)  one  obtains 


Sr(fcli'.Mt')  A(tv)  pev(£',t')  =  p*(x  ,t) 


where 


(30b) 


hi,  (xj  |  xj,  1 1 1')  «  (EXj  (xj  (Xj,t  1 1')  5ji)  , 


(31b) 


nA(2n)  =  A(Ixn)  ;  A(tv)  -(*5  5ji)  ;  aJ«A(tJ)  ;  . 

(33) 

[cf.  equations  (5),  (7),  (1 1),  and  (29).]  The  matrix  operator  A  can  accommodate  a  vector 
dependent  variable  because  it  is,  by  definition,  a  diagonal  matrix.  The  vectorial  form  of  the  elapsed 
time  tv  is  simpler  and,  therefore,  more  convenient  than  its  diagonal  matrix  tv  counterpart  would 
be.  [cf.  equations  (1  lb)  and  (1  lc).J  In  the  wave  approach  the  couplings  are  compatibly  defined  in 
terms  of  a  set  of  junctions.  The  junctions  define  the  insitu  boundaries  and  the  interconnections 
among  the  dynamic  systems.  A  terminal  position  at  a  boundary  (a  junction)  in  the  (j)th  dynamic 
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designated  x»j.  A  terminal  position  vector  may  then  be  defined 


{xtj}  . 


(34) 


The  set  of  the  terminal  position  vectors  define  the  extents  and  the  boundaries  of  the  dynamic 
systems  in  the  complex.  To  complete  the  definition  of  the  model  of  the  complex  it  is  necessary  to 
define  the  junction  matrices 


Jlc  -  (T&)  • 

(35a) 


The  junction  matrix  element  T££ji  defines  the  conversion,  partially  or  totally,  of  an  incident  wave 
propagating  in  path  a,  on  a  junction  at  the  terminal  position  Xbi  in  the  (i)th  dynamic  system,  into  a 
wave  that  emerges  and  propagates,  in  path  v,  away  from  the  junction  at  the  initial  position  xu  in 
the  (j)th  dynamic  system,  see  Figure  4b.  The  conversion  is  ar  ^.^anied  also  by  a  temporal 
integral  operation.  Stating  the  process  just  disc>  ssed  more  explicitly,  with  minimal  loss  in 
generality,  one  may  cast3 


Tvk 


m  m  m 


(35b) 


where  it  is  assumed  that  the  spatial  and  temporal  operative  properties  of  the  junction  matrix  can  be 
vie 

separated;  Abo  is  a  temporal  integral  operator  and  y0*5  is  a  spatial  integral  operator 


3With  certain  modifications  and  extensions  in  the  formalism ,  one  may  account  for  situations  in 
which  the  external  drives  may  be  placed  also  within  the  junctions.  Such  modifications  and 
extensions  were  introduced  in  Reference  2.  However,  such  considerations  lie  beyond  the  scope  of 
the  present  paper. 
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Abo  or  |  t'")  fe  (t'")  =  aS  (0)  f  (t")  ; 

■»  w 

»°‘ (ib  I  x")  =  I  yb(Xb  |x-)d4"  ...  ;  dS"  -  (dij-  S,i)  , 

(36a) 

Vic 

where  Abo  (0)  is  merely  an  algebraic  matrix  factor,  fe  and  f  are  arbitrary,  but  well  behaved, 
functions  of  the  (dummy)  temporal  variables  t"  and  l"'  respectively,  and  yb  and  y*  are  algebraic 

m  m 

matrices  in  the  (dummy)  spatial  vector  variables  x"  and  x"'  and  the  terminal  position  vectors 
Xb  and  Xfc.  To  exhibit  the  properties  of  the  junction  matrix  even  more  explicitly,  the  simple  model 
is  used  yet  again.  For  this  model  the  junction  matrix  is  defined  in  the  form 

Tg  (t£>  -  5k  Ag  80b  , 

■  *  mm  m 

(35c) 


where 


(36b) 


AjSiiCtgii)  =  A^i(0)  ACtgji)  , 

80t«Js(t"-t"')dt'"-  , 


(36c) 


«k-  S  (S'" -  2t)  -  ( 8  (Sj"  -  Sw)  Sj, )  ; 

5°b  -  1  8(xb- 1")  dS". . . .  ;  dx"  *  (dif  8ji)  , 

(36d) 

Aboji  (0)  and  are  merely  algebraic  factors.  As  equation  (36c)  makes  clear,  the  elements  of  Xt£ 
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define  a  set  of  simple  time  delays  that  are  imposed  on  the  appropriate  waves  as  they  traverse  the 
junctions. 


At  this  stage,  it  is  instructive  to  note  in  passing,  a  differential  feature  between  the  initial 
models  for  the  wave  and  the  modal  formalisms.  In  the  wave  approach  the  intial  model  is  defined  in 
terms  of  properly  extrapolated  dynamic  systems;  these  dynamic  systems  are  uncoupled.  A  proper 
extrapolation  is  achieved  by  manipulating  the  model  so  that  die  junction  matrices  vanish;  e.g., 
rendering  the  Abo's  in  equations  (35b)  and  (35c),  identically  zero.  On  the  other  hand,  in  die 

~  J 

modal  approach  the  initial  model  is  defined  in  terms  of  properly  isolated  dynamic  systems;  these 

dynamic  systems  are  uncoupled.  A  proper  isolation  is  achieved  by  manipulating  the  model  so  that 

vie 

the  nondiagonal  elements  in  the  junction  matrices  vanish;  e.g.,  rendering  the  Abo’s,  in  equations 
(35b)  and  (35c),  diagonal  matrices. 


Reiterating  then,  die  model  of  the  complex  is  defined  in  terms  of  die  propagator  matrix  h«  , 
die  terminal  position  vector  x„  and,  finally,  the  junction  matrix  .  Analogous  to  equations  (5) 

and  (1 1),  and  with  the  help  of  Figure  4  and  equations  (30)  and  (32),  one  may  formally  cast  the 
equation  for  the  response  vector  pv  ( x,  t)  that  is  solely  in  path  v  -  -  die  partial  response  vector  -  - 
in  the  matrix  form 


p.(x t)]  -pv<i'.  t) 


(37) 


where  the  partial  impulse  response  matrix  operator  is  of  the  form 1 
6v(x|x',t|t')  -  h*(x|£\t|t')v 

+  2.ZZ  h_(«|£"',t|t'")T£  fi'Cli'.rit')  , 

Jr  h  ^  ^  ^  m  » 

(38a) 


and  the  external  drive  vector  is  of  die  form1 
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(39) 


pe(x,  Pee(x',  t)  ;  C«  (Oj  5ji) 


[cf.  equations  (8)  through  (10).]  With  the  aid  of  equation  (35b),  the  partial  impulse  response 
matrix  operator  may  be  stated  in  the  mote  informative  but  abbreviated  form 


hv  *  h*  v  +  Z  h*  yic  y0b  h° 

-  »  «  kbo  »  «  -  m  m 


(38b) 


It  is  of  interest  to  note  that  equation  (38a  &  b)  is  amenable  to  interpretation  in  a  manner  that  is 
commonly  employed  in  structural  acoustics;  namely,  the  division  of  the  response  into  direct  and 
reverberant  parts.  Of  course,  there  are  a  number  of  different  ways  of  defining  and  cutting  this 
division.  In  establishing  equation  (38a  &  b),  however,  a  division  is  already  suggested.  The  direct 
part  is  to  be  associated  with  that  part  of  the  response  that  occurs  prior  to  interactions  with  any  of 
the  junctions.  Once  the  response  is  contaminated  by  interactions  with  a  junction,  it  is  portioned  to 
the  reverberant  response.  Reverberation  in  this  context  pertains  to  the  initial  and  subsequent 
interactions  of  waves  with  the  junctions  of  the  complex.  [Commonly,  however,  reverberation  is 
reserved  to  situations  in  which  large  multiplicity  of  tins  kind  of  interactions  are  present;  in  fact, 
large  enough  so  that  the  initial  interactions  become  insignificant]  It  is  clear  that  tire  direct  part  of 
the  response  is  associated  with  the  first  term  on  the  right  of  equation  (38a  &  b).  This  term  is 
dubbed  the  direct  term  in  the  impulse  response  matrix  operator.  It  is  equally  clear  that  the  other 
terms  are  associated  with  the  reverberant  part  of  the  response.  These  terms  are  dubbed  the 
reverberant  terms  in  the  impulse  response  matrix  operator.  The  division  is  then  unambiguous  and, 
as  can  be  readily  recognized  and  verified,  remains  so  throughout  the  wave  approach  to  tire 
derivation  of  the  impulse  response  operators.  The  corresponding  equations,  to  equations  (38a)  and 
(38b),  for  the  simple  model  are 
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fiv(8 1  1 1 0  -  Ei(t  |  x',  t  It')  Aft'-)  v 

+  115  EKiil'-.tir)  A(t*)Tg  fi-cx"  1  t"  1  to  , 

(38c) 

hv«  h*  v  +  I  h^A*  80b  h°  ; 

—  mm  kbO  m  m  m  mm 

h* -h*  AV  ;  A^-CA^jiCOJACtbojOS01)  , 

(38d) 

respectively,  [cf.  equations  (35c)  and  (36b)  through  (36d).)  The  notations  are  chosen  so  that  die 
reduction  of  equation  (38a)  to  that  of  equation  (38b)  may  be  instituted  without  confusion;  e.g., 
note  the  cyclic  form  in  with  C  — »  b  — >  k  -4  v;  see  Figure  4b.  In  this  connection,  situations 

may  arise  in  which  the  junction  matrix  may  be  a  priori  of  simpler  form.  For  example,  one 

may  find  that  the  terminal  position  Xb  of  incidence  and  the  terminal  position  x*  of  departure  are 
the  same;  e.g.,  as  depicted  in  Figure  4c.  Then  Tj£  -4  T$  6*b  ;  •  Tj^.  The  summation  in 

equation  (38),  and  in  subsequent  equations,  may  then  be  accordingly  and  a  priori  condensed.  It 
may  also  be  noted  that  the  summation  over  die  terminal  positions;  e.g.,  in  equation  (38),  implies 
discreteness.  If  the  discreteness  is  found  to  be  inappropriate,  the  summation  may,  in  part  or  in 
whole,  be  replaced  by  the  appropriate  integration.2  It  may  be  proposed,  at  this  stage,  that  the 
following  notations  be  recorded  for  subsequent  use 

A/b«  A(T?b)  «  (AC&pSjj)  ;  t)bj  -  X^xjlXbj)  ; 

Av«/«  A(t*,)  -  (A(x^j)6ji)  ;  x^-j«  xj(xlj|xj0  ; 

A«b*  A(x  tb)  *  (A(Xtbj)8ji)  ;  X ibj  *  Xj  (Xgj  |  Xbj)  . 

(40) 

Equation  (38)  appears  impure  in  the  sense  that  a  number  of  h°'s  show  in  die  reverberant 
terms  of  the  expression  for  h.  Means  need  be  devised  to  purify  equation  (38).  For  this  purpose 
equation  (38b)  is  manipulated  to  construct  the  relationship 
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where 


(41a) 


(42a) 


-  y0b  hZ,  y  f  A?$  ;  S;e  -  (6;o  fy)  , 

*  ik  ™  ib  *  * 


(43a) 


cvk 

and,  again,  note  the  cyclic  form  in  Abo  with  c— »b— »k— »v—»c;see  Figure  4.  Substituting 
equation  (41a)  in  equation  (38b)  yields 

hv  *  h*  v  +  (  £  h*  y*  A^)  Do  {y0b  h°  o 

*  »  «  kbo  *  m  ■»  Li  «  «• 


+  ( S  A^f)  kj,  5}  . 


(44a) 


where 


D$«(BS)-1  . 


(45a) 


The  reverberant  terms  in  equation  (44a)  are  clearly  pure,  and,  therefore,  so  is  £v.  Thus,  the 
construction  of  die  Bo's ,  stated  in  equation  (41a),  achieved  their  designed  purpose.  The 
construction  of  the  matrix  Bo  is  central  to  the  reverberant  part  of  the  response.  Whenever  a  set  of 
parameters  in  this  matrix  renders  (the  magnitude  of)  its  eigenlike  value  small  compared  with  unity, 
die  reverberant  response  is  high.  A  highly  reverberant  response  is  potentially  indicated  if  die 
eigenlike  values  of  Bo  are  small  and  numerous.  Indeed  the  potential  necessity  for  die  existence  of 
such  small  eigenlike  values  in  B§ ,  prompted  its  construction.  In  die  frequency  domain,  a  set  of 
parameters  that  yield  high  reverberant  response,  as  just  prescribed,  is  dubbed  a  set  of  ’'resonance" 
parameters.  A  common  example  of  a  parameter  in  this  set  is  the  resonance  frequency.  It  is 
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expected  that  invariably  the  contributing  terms  to  die  reverberation  in  die  response  will  harbor  die 
matrix  factors  D&  that  are  the  inverse  of  theD&'s.  [cf.  equations  (44a)  and  (45a).]  It  is  noted  that 

m  m 

of  bof 

the  matrix  B&  is  a  temporal  operator  as  Ae£  or,  equivalently,  Ae$  is.  In  die  spatial  domain  this 
operator  is  neuter.  This  statement  is  clarified  noting  that  y 1 °b  h*  y  f  is  not  a  spatial  matrix  operator. 

m  m  m 

Since  B§  is  a  temporal  matrix  operator  so  is,  of  course,  D& .  However,  whereas  B&  is  a  proper 
operator,  its  inverse  may  need  special  attention  when  manipulated.  Nonetheless,  the 
expression  for  the  impulse  response  matrix  is  pure  in  the  sense  that  it  is  a  functional  purely  and 

solely  of  the  properties  of  die  dynamic  systems  and  the  couplings  among  them;  namely,  of  the 
h«  s ,  xt' s,  and  ’ s .  Again,  making  use  of  the  simple  model  yields,  for  equations  (41a)  through 
(43a),  the  corresponding  equations 


(41b) 


B^  -=  [  I  -  )  ( .X  Ac;  )  X  5cbX  6;c]  , 

-  80b  h-  5f  ;  8;o  -  (8;0  5ji)  , 


(42b) 

(43b) 


respectively.  Substituting  equation  (41b)  in  equation  (38c)  one  obtains  for  the  simple  model 

Sv-h*  v  +  (X  hi  Bt  aS)  nS  {60b  hi  o 

n  ■,  m  kbC  •  «  ■>  »  m 


*(Z  a5?)  80*  *1  5}  ;  hi-E''Av  , 

fet  *  -  -  «  -  - 


(44b) 


where 


pS-OJ&r1 


(45b) 


As  in  equation  (45a),  it  is  noted  that  the  matrix  B&  in  equation  (45b),  is  a  temporal  operator  as  Ae$ 
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is.  Lflce  its  counterpart  in  equation  (45a),  in  the  spatial  domain  this  operator  is  neuter.  This 
statement  is  clarified  noting  that  80b  h£,  8f  is  obviously  not  a  spatial  matrix  operator. 

Once  equation  (44)  is  derived,  die  grand  equation  for  die  response  vector  p(x,  t)  of  the 
complex  may  be  stated  in  terms  of  die  impulse  response  matrix1 

h(£|£',t|t')-Xhv(£|i',t|t')  ; 

h(  1 14',  t  It')  -  X  X  {hX(i  |x',  t  It') v  S„,  +  hS(x  |t'.  t  It')  o}  . 

*  ^  VO  *  at  ai  •  ■ 

(46) 


and  the  external  drive  vector1 


(39) 


in  the  form 


(47) 


[cf.  equations  (8)  through  (10).]  The  expression  for  the  impulse  response  matrix  operator  h  in 

A 

die  second  of  equation  (46)  is  of  some  interest  Hie  double  summation  in  h  is  introduced  to 
accommodate  die  detail  form  of  die  drive  Pe  *  ^  The  first  term  in  the  curly  brackets  is  the 
direct  term  and  the  second  term  describes  the  reverberant  terms  in  the  impulse  response  matrix 
operator.  The  direct  term  in  the  impulse  response  matrix  operator  is  always  readily  distinguishable 
from  die  reverberant  terms;  equation  (46)  is  a  prime  example  of  this  statement  [cf.  equation  (38).] 


Equation  (47)  is  expressed  in  the  temporal  domain  and  is  suitable,  if  not  essential,  to 
ascertaining  die  response  of  a  complex  to  transient  external  drive  systems  -  -  external  drives  that  are 
localized  in  the  temporal  domain.  Inserting  equation  (44)  in  equation  (47)  results,  even  in  models 
of  simple  structural  complexes,  in  an  equation  that  is  cumbersome  and  difficult  to  interpret  Some 
relief  may  be  derived  if  the  complexes  are  assumed  to  be  stationary  in  the  temporal  domain.  This 
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kind  of  stadonarity  is  imposed  on  the  complex  by  assuming  that  die  dynamic  systems  and  the 
couplings  among  them  remain  unchanged  with  the  passage  of  time;  i.e., 

h'  (i|i\  t  It') ■  ait)-1'2  ; 

tm  ~  ~  *  “ 

aS  (no  *  a tt)-i/2  a2  (t^—  t"')  . 

(48a) 

[cf.  equations  (36a)  and  (38a).]  In  the  simple  model  the  conditions  of  stationarity  are  more 
stringently  imposed;  namely, 

h*  (x|x',t|t')»  (2JC)-1'2  h*  (£' |  x,  t - t')  A  (xv)  ; 

/£  -  «,  (0)  A  d&ji)  S“)  , 

(48b) 

and  the  A^ji  (0)'s  and  the  s  are  merely  algebraic  factors,  [cf.  equations  (36b)  and  (38c).] 
Substituting  equation  (48)  in  equations  (42),  (44),  and  then  in  equation  (46),  helps,  to  be  sure,  but 
die  formalism  still  remains  cumbersome  and  difficult  Further  relief  may  be  derived  if,  in  addition, 
the  equation  far  the  response  is  Fourier  transformed  from  die  temporal  domain  into  the  frequency 
domain.  The  equation  in  die  Fourier  domain  is  usually  more  suitable  for  ascertaining  die  response 
for  "steady  state"  drives  and  for  complexes  that  are  stationary  with  respect  to  the  temporal  domain. 


IMPULSE  RESPONSE  OPERATORS  IN  THE  FREQUENCY  DOMAIN 

To  this  point  die  analytical  procedures  and  developments  were  conducted  in  the 
{x ,  t } -space.  Is  there  any  advantage  in  transforming  the  equations  derived  in  die  preceding 
section  into  other  related  spaces;  e.g.,  by  a  Fourier  transformation?  It  was  argued  in  die 
introduction  that  such  a  transformation,  into  a  specified  domain,  may  be  advantageous  when  die 
complex  is  stationary  with  respect  to  the  variables  that  define  this  domain.  In  this  section,  only  a 
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temporal  stationarity  is  considered.  If  the  dynamic  systems  and  the  couplings  among  diem  remain 
unchanged  with  die  passage  of  time,  as  and  stated  in  equation  (48),  then  under  these 
crtnAinrtnn  equation  (30)  may  be  recast,  by  a  Fourier  transformation,  in  the  frequency  domain 


hI(x|{',(o)p„(5',«i) 

m  ** 


pI<X,  CD  ) 


(49a) 


or  in  the  abbreviated  form 

hi  (to)  Pev  (co)  *  pi  (co)  , 

m  ~~ 


(49b) 


where 

hl(x|x',©)  ->  (2jc)*1/2  hi (it |x',  t  - 1' )  exp[-ico(t-t')]  , 

p«r(x,a>)  -♦  Jdt'peo(x,t')exp(-i©t')  , 

pl(x.co)  -»  (2ti)-1/2  J dt  pl(x,  t)  exp(-icot)  . 

[cf.  equation  (1).]  It  is  also  noted  that  the  arrows,  replacing  die  equality  signs,  indicate  that  die 
quantities  are  related  by  Fourier  transformations.  However,  unlike  the  notations!  procedure 
adapted  in  equations  (1)  through  (3),  the  letterings  remain  the  same.  This,  together  with  the 
reduction  of  equation  (49a)  to  equation  (49b),  is  instituted  to  keep  the  notations,  which  are  already 
unwieldly,  to  a  mfnimnm.  The  relief  in  transforming  from  the  time  domain  into  the  frequency 
domain  is  achieved  mainly  because  the  temporal  operators  in  the  temporally  stationary  formalism 
are  eigen-operators  with  respect  to  the  Fourier  eigenfunction  (2k)*  1/2  exp  (icot).  Therefore,  in  the 
frequency  domain  these  temporal  operators  are  simply  replaced  by  their  eigen-values.  These  eigen¬ 
values  are  merely  algebraic  factors.  This  replacement  is  particularly  explicit  in  the  formalism  of  the 
simple  model  For  the  simple  model  equation  (50a)  can  be  cast  in  die  form 


(50a) 

(51a) 

(51b) 
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hl(x|£',a»  -» h£,(£|x',t-t')  A  (tv)  (2*)'ly2  exp[-ia>(t-t')]  ; 
hH(x|x',a>)  «Bh*(£|$',a>)V(©;cv)  mhl«o)  ; 

«■  m  m  m 

A(tjj)  V(O)Tjj)  ) 

[(2jc)*1/2  exp  (icot)]  =  J  [(2n)*1/2  exp  (icot)]  ; 

F[A(tJ)]  F[V(tDtjj)]J 

V(tDTv)  =  (expC-icotjj  )5ji)  =  (V(cotjj)8ji)  . 

(50b) 

where  F  is  an  arbitrary,  but  well  behaved,  functional  of  the  dependent  quantity.  Imposing  the 
temporal  stationarity,  stated  in  equation  (48),  on  equations  (42)  and  (43),  one  obtains 

Bo  (to)  =  [  I  -  (X  ASf(co))  (X  A*f(e>))X  8cb  X  5ra  ]  , 

^  4  (52) 


where  for  the  more  general  me  ■ud 


Ajf(©)  -  h^f(to)  AS(CO)  ; 

A$(cd)  -►  (2jc)*1/2  A$(t"-0  exp[-i(D(t"-t"')]  ; 
h^bf  -  y0b  h®(to)  yf  , 

*  jk  m  m 

(53a) 
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and  for  the  simple  model 


Ajf(co)  =  h£bf<«)  aS(o»  ; 

A$«d)  *  ( A^|ji  (0)  V(cox^ji))  ;  Vtou^j)  =  expC-icoi^i)  ; 
h2w(«o)  =  h2*f(o>)  V(«t?f)  ;  h^bf  (co)  -  80b  E®(co)  8f  ; 

V(tot^f)  =  (exp(-ico^fj)  8ji)  *  (V  (artg^)  Sji)  . 

(53b) 

[cf.  equation  (40).]  As  expected  and  anticipated,  equation  (53b)  is  more  explicit  and  simpler  than 
equation  (53a).  Imposing  die  temporal  stadonarity,  stated  in  equation  (48),  on  equations  (44a)  and 
(44b)  yields 

fiv(co)«=  h*  (cd)  v  +  X  hX(<o)  yk  a£(g»  D&(co) 

m  mm  KDO  m  m  m  m 

{y0b  h®  (cd)  o  +  X  A^f  (CD)  yOe  h^(CD)  5  }  , 

m  m  m  fe£  “  m  m  m 

fiv(CD)  «  h^(co)  v+  X  hX(co)  8k  a£(cd)  D&(cd) 

m  mm  kbO  »  m  m  m 

{80b  h®  (cd)  a  +  X  A$f  (co)  8°*  h^(a»  5}  , 

-  -  -  i*  •  •  “ 

respectively,  where 

D&(co)  -  [B^cd)]-1  • 


(54a) 


(54b) 


(55) 


[cf.  equation  (45).]  It  is  noted  that  the  matrix  B&  (cd),  unlike  B&  in  equation  (42),  is  merely  an 
algebraic  matrix  factor,  and  therefore  so  is  Do  (co).  Consequently,  in  the  frequency  domain,  the 


impulse  response  matrix  h  (to)  is  of  simpler  structure  than  it  is  in  the  temporal  domain.  Again,  the 
direct  term  hi  v  in  the  impulse  response  matrix  operator,  in  equation  (S3  a&b),  is  clearly 
distinguishable  from  the  reverberant  terms.  Moreover,  it  emerges  that  the  grand  equation  for  the 
response  vector  p  (x,  <o )  of  a  model  that  is  temporally  stationary  is,  as  expected,  algebraic  in  die 
frequency  domain 


[fi(x|x',0))  pe(x',  CO)]  =  p(x,  CD) 


(56) 


where 


h(x|x',co)  *  X  hv(x|x',to)  ; 

sat  **  ^  v  as  ^  ^ 

fi(x|x',  to)  {hl,(x(x/,  co)v  5vo+  ho  ( x  |  xr,  co)  a}  , 

m  ~  ~  V  G  mb  ^  sb  mt  ■  ^  ^  m 

(57) 


Pe(x',tO)  cpeo(t',0))  , 


(58) 


are  the  impulse  response  matrix  and  the  external  drive  vector,  respectively,  expressed  in  die 
frequency  domain,  [cf.  equations  (39),  (46),  and  (47).] 

The  relief,  expected  from  the  stationarity  of  die  complex  in  the  temporal  domain  and  die 
Fourier  transformation  of  equation  (47),  is  materialized  in  equation  (56).  Nonetheless,  as  indicated 
by  equations  (53)  and  (54),  substantial  cumbersomeness  and  difficulty  still  remain  in  equation  (56). 
This  stems,  as  it  does  in  equation  (47),  from  die  multi-spatial  dimensionality  of  the  complex.  How 
much  relief  may  then  be  reasonably  gained  by  restricting  the  model  of  die  complex  to  a  composition 
of  coupled  one-dimensional  dynamic  systems? 
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COUPLED  ONE-DIMENSIONAL  DYNAMIC  SYSTEMS 


As  argued  in  the  introduction,  if  an  impulse  response  operator  is  stationary  in  a  given 
dependent  variable,  it  can  be  rendered,  by  a  Fourier  transformation,  algebraic  in  this  variable. 
Then,  by  focusing  on  a  single  value  of  die  Fourier  conjugate  of  the  variable,  the  equation  for  the 
response  is  effectively  reduced  in  the  dimensionality  of  this  variable  [3,4].  Thus,  if  the  complex  is 
stationary  with  respect  to  £2  in  the  spatial  domain  x  =  {x,  X2},  the  equation  for  the  response  can 
be  reduced,  as  just  specified,  to  account  for  a  model  of  a  complex  composed  of  spatially  one¬ 
dimensional  dynamic  systems,  where  conveniently  xi  s  x,  and  the  reduced  equation  is  specified  at 
a  specific  value  of  k2;  k2  being  the  Fourier  conjugate  of  £2 .  In  this  paper,  for  the  sake  of  brevity, 
this  one-dimensionality  of  the  dynamic  systems  is  imposed  directly  on  equations  (42)  through  (44) 
and  on  equations  (51)  through  (53).  The  reference  to  the  other  spatial  dimensionalities,  if  present, 
is  simply  suppressed  and  omitted  in  the  Fourier  transform.  The  authors  previously  presented,  in 
References  4  and  5,  an  analogeous  analysis  of  the  response  of  a  model  composed  of  multiple 
coupled  one-dimensional  dynamic  systems.  In  these  previous  considerations,  however,  the  one- 
dimensionality  of  the  dynamic  systems  is  imposed  a  priori.  In  the  first  reference,  the  analysis  is 
conducted  a  priori  in  the  {x,  co}-spacc.  In  the  second,  and  the  more  recent  reference,  the  analysis 
is  conducted  a  priori  in  the  (x,  t}-space.  In  this  way  the  relationship  between  subsequent 
considerations,  in  this  section,  and  the  previous  considerations,  in  References  4  and  5,  can  be 
established  with  the  insights  that  such  a  comparison  may  provide  to  the  modeling  and  the 
formalism. 

A  complex  consisting  of  one-dimensional  dynamic  systems  is  sketched  in  Figure  5.  The 
complex  is  appropriately  defined  in  terms  of  just  two  junctions,  junction  r  and  junction  q  [1-5]. 
Thus,  the  model  is  specified  in  terms  of  two  propagator  matrices,  two  terminal  vectors,  and  two 
junction  matrices.  For  the  more  general  model,  die  propagator  matrices  are  in  the  forms 
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h“(x|x',t  |t')  =  (hS,  (xj  | xj,  t  |t')8ji)  ;  x  «  {xj}  ; 

x,  -»  x  ;  xij  Xj  ;  a  *  r  or  q  , 

‘-1  (59a) 


and  for  the  simple  model  in  the  forms 

h“(x|x\t|t')  *  h“(x|x',t|t')  A(t°)  ; 

m  ~  ~  at  ^  ^  «  ** 

j£<x|x',t|t')  «  (h“(xj|xj ,  tlOSji)  ; 
lo(x,|x,0  =  {'tjI(Xj|xp}  =  {i§}  . 


(59b) 


The  two  terminal  position  vectors  are  specified  in  the  forms 

Xo  »{xaj>  ;  a  =  r  or  q  , 

These  terminal  position  vectors  are  common  to  the  mare  general  and  the  simple  models.  The  two 
junction  matrices  far  die  more  general  model  assume  die  forms 

T£  =  5bk  Sbo  S«w  (f  —  Svo)  T0  ; 

tn  IB  Ms  Ml  Ml  Ml 

To  =  6  (x  "  -  Xo)  Ac  J  5  (Xo  -  x  "')  d  x'"- . .  ; 

**  Ml  ~  *  Ml  ~  Mi 

A o  *  (Aoji)  ;  a  =rorq  ; 

Ao(t"lOfe(t/w)-A«(0)f  (t")  , 

(61a) 


and  for  the  simple  model  they  assume  the  forms 
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«  8bk  Sbo  8oa  (I  -  Svo)  Ta  ; 

m  m  m  m  m  m  m 

Ta  =5  (x -xa)Aa  |  5  (Xa-x)dx'...  ; 

Aa  =  (Aaji  (Xaji)  8  )  J  Aaji  (Xaji)  *  A«ji  (0)  A  (Xaji)  . 

(61b) 

[cf.  equations  (30)  through  (36).]  Substituting  equation  (59)  and  (61)  in  equations  (42),  (46),  and 
(44)  yields 

Ba=[i-Af  4]  ;  Da*  (Ba r1  ; 

w  *= r  as  as  as 

4-8°%!,S(,A„  , 

m  m  *:  m  m 

(62) 

fi  fia  ;  fia  =  (h“+hg)  a  +  h|  p  , 

m  I,q  a,  m  m  m  ms  mV 

(63) 

hg  =  h“  8p  ApDp  A|5°°h“  ; 
hJ-hiSpApDp^hi  , 

(64) 

0 

respectively.  It  is  noted  that  the  matrix  quantities  Ac  and  A«  are  temporal  operators  and,  therefore, 
so  is  Ba  and,  consequently,  so  is  Da.  In  this  sense  equations  (62)  and  (64)  cannot  be  manipulated 
cavalierly;  commutations  and  expansion  procedures  must  be  carefully  viewed  and  performed. 
Equations  (59),  (61),  and  (62)  through  (64)  are  die  equations  for  die  response  of  die  model  of  a 
complex  composed  of  one-dimensional  dynamic  systems.  Compared  with  equations  (35),  (36), 
and  (41)  through  (44)  for  die  multi-dimensional  dynamic  systems,  the  reduction  achieved  in 
equations  (59),  and  (61)  through  (64)  is  dramatic.  Moreover,  it  is  noted  that  the  spatial  operative 
properties  of  the  junction  matrices  in  equations  (59a)  and  (61a)  and  those  in  equations  (59b)  and 
(61b)  are  identical;  a  further  and  considerable  simplification  indeed.  The  difference  in  the  junction 
matrices  described  in  these  two  sets  of  equations  lies  entirely  in  die  temporal  domain,  but  the  one  in 
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equation  (59b)  and  (61b)  are  by  far  the  simpler.  Finally,  it  is  observed,  from  equations  (63)  and 
(64),  that  the  impulse  response  matrix  operator  ,  for  the  model  composed  of  coupled  one¬ 
dimensional  dynamic  systems,  harbors  six  (6)  distinct  terms:  two  (2)  direct  terms; 
h“  a,  a  =  r  and  q,  and  four  (4)  reverberant  terms;  h§  oc  and  h£  f$,  a  =  r  and  q.  Although  the 
impulse  response  matrix  operator  for  the  one-dimensional  model  is  significantly  simpler,  some 
further  simplification  in  die  model  is  still  required  to  attain  reasonable  descriptions.  The  first  in  this 
process  is  to  bring  to  bear  the  temporal  stationarity  onto  the  one-dimensional  model.  Thus, 
imposing  the  stationarity  in  the  temporal  domain  on  the  one-dimensional  model  recasts  equation 
(48a)  in  the  form 

h“  (x|x',t  |t')s  (2tc)-1'2  h“  (x  |  x ',  t  -O  ; 

m  m  ~  ~ 

Aa  (t"|t'")  s  Ac  (t"~  t'")  , 

(65a) 


and  for  the  simple  model,  equation  (48b)  yields 

h“(x|x',t|t')  s  (2it)”1/2  h“(x'|x,t-t')  A(t“)  ; 

Ac  =  (Aaji(O)  A(Xoji)  6  )  » 

(65b) 

and  the  Aa  (0)'s  and  the  xa's  are  algebraic  factors.  Substituting  equations  (65a)  and  (65b)  in 

m  m 

equations  (62)  through  (64)  yields  further  simplification.  However,  the  significant  simplification 
is  attained  when  the  resulting  equations  are  transformed  into  die  frequency  domain.  Such  Fourier 
transformation  yields 

Ba«a)  «  [I-Ap«o)  a£(o»]  ;  Da(Q»  *  [Bafco)]*1  ; 

Aa(CO)  =  8°P  hl(tO)  8a  Aa(CO)  «  gL(£pl*a,C0)  Ao(CO)  , 

(66) 
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£(®)  “  #5  *  fia(«)  *  [h~  (CO)  +  hg  (CD)]  a  +  hp(co)£  , 

hg(co)  -  h“«o)  8p  Ap(co)  Dp(co)  A^(co)  8°°  h“ (a» ; 
hp(a>)  =  h“(a»  8p  Ap(co)  Dp(co)  8°^  l£(a»  , 

respectively,  where  for  the  general  model 

gpQtpl£a,CD)-4  (In )~m  80P  hp(x  |x',  t-t')  8a  exp  [-im(t-t')] 

m  m  m  m 

Aa(CO)  ->  (27C)~ 1/2  Aa(t"-  t'")  exp  [  -i<M(t"-  t'")]  . 

*1  m 


(67) 


(68) 


(69a) 


and  far  the  simple  model 

gLQcpl£a,CD)  «  g!(Xp|Xa,CO)  V(0)tP  )  ; 

•  M  *  ’ 

sLCiplio.®)  -*  Sw  6aMp[-iO)(t-t')]  ; 

wm  mt  mt  mt 

Y  (cox^)  «  (exp  [  — ico  xf(xpi  |  x^)]  Sp)  ; 

Aa(<o)=  (Aaji(0)V(anaji))  ;  VCcoiaji)  =  exp(-imtoji)  . 

(69b) 

where  g* (co)  and  g“(co)  are  die  impulse  response  functions  that  are  associated  with  die  impulse 

m  m 

response  operators  h“(co)  and  hi(03),  respectively,  [cf.  equations  (4)  and  (40).]  Comparing 
equations  (66)  through  (69)  with  equations  (50)  through  (54),  again,  reveals  the  considerable 
simplification  attained  in  modeling  the  complex  in  terms  of  coupled  one-dimensional  dynamic 
systems.  However,  to  make  ready  contact  with  the  formalisms  developed  and  used  in  Reference  5 
and  References  1  and  2,  further  simplifying  assumptions  need  to  be  imposed  on  the  one- 
dimensional  model  of  the  complex. 
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BASIC  COUPLED  ONE-DIMENSIONAL  DYNAMIC  SYSTEMS 

Using  equations  (65b)  and  assuming  the  coupled  one-dimensional  dynamic  systems  to  be 
basic,  one  may  explicitly  define  and  impose 

(2jc)"1/2  h“(x|x',t-t')«ga(x|x')  50t(t-t')  ; 

Aa+(xa)  *  U^Cx-xO  A(t°)  ;  ga(x|x')  e  U°+(x-x')  ; 

**  m  m  m 

t-  {(Cfr1 1  xj  -  xj  i }  ;  cf -qti-in?)-1  ; 

-  xj)  -  U[(Xj  -  xj)  sign  (Xoj  -  XK)]  ; 
a-rorq  ;  £«=rorq  ;  a*p  ,  (70) 

where  Cqj  is  a  temporally  independent  speed  of  propagation  toward  junction  a  and  Tlf  is  the  loss 
factor  associated  with  this  propagation  in  the  (j)th  dynamic  system  [3].  [cf.  equations  (6b)  and 
(59).]  Consulting  equation  (41),  and  using  equations  (62),  (65b),  and  (70),  one  obtains 

Ba.  [I-a|4]  ;  p0-(Ba)-’  ;  4°*($).A(&)  ; 

(71) 

Moreover,  from  equations  (63),  (64),  (65b),  and  (70)  one  obtains 

h°  «  J  {[A^t*)  +  g£*(x |x')]  o  +  g?p(x |x')  p}  dx'  80l(t-t') ...  , 

m  m  m  mm  m  m 

(72) 

*5,(51*')  -  A°*<i5)  *t  Dp  A*^)  A„  A^iS,)  ; 
jifeli')-  A°*(j£>  Ap  Dp  A*(£)  . 

(73) 
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where  Dp,  in  equation  (73),  is  as  stated  in  equation  (71),  and  ?£/,  xjp,  and  x£p,  are  deduced  from 
equation  (40).  Equations  (39),  (47),  (63),  and  (71)  through  (73)  are  employed  in  Reference  5  to 
exemplify  the  nature  of  the  response  vector  p  (x ,  t)  of  a  model  of  a  complex  composed  of  multiple 

mt 

coupled  basic  one-dimensional  dynamic  systems.  In  this  reference,  however,  die  one- 
dimensionality  of  the  basic  dynamic  systems  was  assumed  a  priori. 

Using  equations  (69b)  and  assuming  the  coupled  one-dimensional  dynamic  systems  to  be 
basic,  one  may  explicitly  define 

h°(x|x',0»  «h°(x|x')  ;  V°+(coTa) -  U°+(x-x')  V(©t°)  ; 

Nt  m  mt  Mr  mt  m  Mr  Mr  Mr 

T°  *  {( Cf) 1  |xj-  Xj  |}  . 

[cf.  equation  (6b).]  From  equations  (66)  and  (69b )  one  obtains 

B«(©)  -  [I-Ap(m)  A £(©)]  ;  Da(Cf>) « [Ba((n)]'1  ; 

Aa(0»  -  ^(a>)  ;  Ipa  *  (cxpHfiO^)^)  -  (V(‘4xj)5ji)  • 

(75) 

[cf.  equations  (40)  and  (71).]  Moreover,  from  equations  (67),  (68),  (69b),  and  (74)  one  obtains 
ha(o>)«J  {[V^cot0)  +  g^(x|^©)]a  +  g®  (x|x',a>)P}  dx'--.  , 

^(xfx'.O)  -  Vjp  Ap(£0)  Dp  (a)  vj.  Aa(»)  ; 

(77) 


where 
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r  *  Y0+(“Ia)  ■  (U?+  (xj  -  xj)  exP  (-  icoxij)6ji)  ; 

V/J+  *  ^(cwj)  =  (Up  (xj  -  xpj)  exp  (-  icot/pj)6ji)  ; 

VS/+  =  Va+(a>t£/)  *  (U?+  (xaj  -  xj)  exp  (-  icoTa/j)8ji)  ; 
v£p  -  =  (U?+  (x^  -  xpj)  exp  (-  ionapj)8ji)  . 

(78) 

[cf.  equation  (40).]  Equations  (56),  (58),  (67),  and  (75)  through  (77)  are  identical,  except  for 
minor  notational  differences,  with  the  equations  derived,  in  References  2  and  4,  for  the  impulse 
response  matrix  operators  and/or  functions.  In  these  references,  however,  the  one-dimensionality 
of  the  basic  dynamic  systems  was  assumed  a  priori.  It  may  be  pointed  out  that  in  the  frequency 
domain,  temporally  stationary  models  of  complexes  composed  of  one-dimensional  dynamic 
systems  can  be  treated  with  relative  ease,  even  if  these  dynamic  systems  are  not  basic.  Thus,  for 
example  C*\  in  equation  (74),  may  be  assumed  to  be  frequency  dependent,  rendering,  in  turn, 
frequency  dependent  Also,  more  complex  propagation,  such  as  may  occur  in  a  model  that  befits 
cylindrical  geometry,  may  be  treated  in  the  frequency  domain  more  readily  than  in  the  temporal 
domain  [4].  However,  in  these  instances  one  needs  to  return  to  the  more  general  equations;  the 
equations  for  the  simple  model  cannot  quite  handle  such  situations.  Nonetheless,  the  formalism 
herein  stands  ready,  with  some  increase  in  cumbersomeness,  to  accommodate  even  these  deviant 
situations.  This  increase  in  cumbersomeness  is  miniscule  compared  with  that  introduced  by 
increase  in  dimensionality.  As  is  clearly  demonstrated  in  the  wave  approach  just  discussed: 
increase  die  dimensionality  and  the  task  is  increased  by  much  more  than  a  linear  proportionality. 
This  leads  to  a  few  final  remarks.  The  modal  formalism  appears  less  dependent  on  the  spatial 
dimensionality  of  the  dynamic  systems;  see  Section  IV.  Of  course  this  is  misleading.  The  spatial 
dependence  of  the  formalism  is  removed  and  stored  away  in  terms  of  sets  of  eigenfunctions. 
However,  the  eigenfunctions  are  strongly  influenced  by  the  spatial  dimensionality  of  the  dynamic 
systems.  Indeed,  making  die  eigenfunctions  for  all  the  dynamic  systems  available  for  storage  is 


43 


the  greater  part  of  the  solution  to  the  spatial  dependence  of  die  response  of  these  dynamic  systems. 
The  availability  and  the  storage  of  die  eigenfunctions  introduce  severe  limitations  on  die  versatility 
of  the  modal  formalism.  An  eigenfunction  is  usually  fixed  in  the  format  of  "standing  waves" 
which  gives  it  quasi-symmetric  properties  with  respect  to  the  spatial  dependence.  If  the  response 
does  not  follow  spatial  patterns  that  can  be  described  by  a  few  of  die  eigenfunctions,  the 
cumbersomeness  of  the  description  may  become  forbidding.  Indeed,  for  a  lopsided  response, 
such  as  may  be  caused  by  excessive  damping,  significant  disparity  in  the  junctions,  substantial 
rectification  in  the  propagation  in  some  of  the  dynamic  systems,  etc.,  the  modal  formalism  may  be 
overtaxed.  Far  example,  try  and  describe  the  direct  part  of  the  response  with  the  modal  formalism! 
On  die  other  hand,  the  wave  formalism  does  not  suffer  this  "locked-in"  limitation;  e.g.,  the  direct 
part  of  the  response  is  readily  determined!  However,  die  spatial  dependence  in  the  wave  formalism 
is  built  piece-by-piece.  Therefore,  the  formalism  is  cluttered  with  die  spatial  dependences.  A 
reduction  in  the  spatial  dimensionality  of  the  dynamic  systems  is  communicated  directly  and 
dramatically  to  the  wave  formalism.  Thus,  although  the  two  formalisms  are  formally 
complimentary,  they  may  also  be  supplementary.  That  is,  in  addition  to  furnishing  the  same 
information,  they  may  also  reveal  different,  yet  useful,  information  in  trying  to  resolve  a  problem 
in  which  the  issue  is  the  response  of  an  externally  driven  complex. 
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Fig.  1.  Schematics  of  a  complex  on  which  are  shown  two  typical  paths;  path  v  and  path  a, 
from  the  localized  external  drive  positions  at  {  x',  t'}  to  the  localized  observation 
position  at  {  x,  t}. 


Fig.  2.  Schematics  of  a  complex  subdivided  into  (three)  dynamic  systems.  The  gj  is  the  self- 
impulse  response  function  of  the  (j)th  dynamic  system  and  gjt  is  the  coupling  impulse 
response  function  between  the  (i)  and  the  (j)th  dynamic  systems. 
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Fig.  3.  Schematics  of  a  complex  modeled  by  coupled  (three)  dynamic  systems.  The  Zj  is  the  self 
impedance  operator  of  the  0th  dynamic  system  and  z^is  the  coupling  impedance  operator 
between  the  (i)th  and  the  0th  dynamic  systems.  The  self  impedance  operator  Zj  defines  a 
set  of  eigenfunctions  ^  and  a  set  of  corresponding  eigen-values  z^. 
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Fig.  4.  Schematics  of  a  complex  subdivided  into  (three)  dynamic  systems.  Propagation  is 

initiated  in  path  a  at  the  position  xj  in  the  (i)th  dynamic  system.  The  propagation  via  this 
path  reaches  the  boundary  at  position  Xbi*  A  component  of  the  propagation  enter  the  (j)th 
dynamic  system  at  the  position  x^j  and  propagates  in  path  n  to  the  observation  position 
Xj  in  the  (j)th  dynamic  system.  [The  latter  propagation  may  continue  until  it  reaches  and 
is  incident  on  the  position  xCj.  ] 


Fig.  4b.  Incident  position  xu  and  exit  position  x^j  are  not  coincident 


Fig.  4c.  Incident  position  Xbi  and  exit  position  xy  are  coincident 
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Fig.  5.  Schematics  of  a  complex  modeled  in  terms  of  multiple  one -dimensional  dynamic  systems. 
Two  junctions  are  necessary;  junction  r  and  q.  These  junctions  are  defined  by  the  terminal 
positions  xa  and  by  the  junction  matrices  A«;  a  *  r  and  a  *  q. 
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